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© Parallel pseudo-random generator for emulating a serial pseudo-random generator and method for 
carrying out same. 

© A parallel pseudo-random generator for emulating a serial pseudo-random generator that generates serial 
outputs such that the next serial output value is based upon an Exclusive OR combination of at least two 
preceding serial output values the maximum preceding serial output value defined as the Pth preceeding serial 
output value, where P is an integer greater than one; comprising: 

A) at least P latches, each latch having an output having a logic value 1 or 0 and an input operable 
upon receipt of a clock signal, for receipt of data for controlling the next logic value on the latch output; 

B) at least P Exclusive OR gates, each having at least two inputs and one output, each Exclusive OR 
gate output connected to a corresponding input of one latch so as to define the next value of the latch 
output upon receipt of the next clock signal; and 

C) means for connecting each input of each Exclusive OR gate to one latch output so that the output 
of each Exclusive OR gate represents the corresponding next value of the latch to which This Exclusive Or 
gate output is connected. 
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PARALLEL PSEUDO-RANDOM GENERATOR FOR EMULATING A SERIAL PSEUDO-RANDOM GENERATOR 

AND METHOD FOR CARRYING OUT SAME 



TECHNICAL FIELD 



The present invention relates to a circuit and associated method for emulating the output of a serial 
5 pseudo-random generator (PRG) or scrambler by a parallel implementation comprising a plurality of outputs 
which represent successive serial outputs of the serial PRG. The invention has particular use in tele- 
communications, where high speed data streams are combined with a serial PRG so as to insure proper 
clocking and for potential security of the data stream. Due to the high-speed nature of such telecommunica- 
tion data, serial PRG's cannot be implemented using complimentary metal oxide silicon (CMOS) circuitry. 
70 Thus there is a need for emulating the serial PRG so that the clock rate of the circuitry is within the 
operating frequency of CMOS circuitry. 



BACKGROUND OF THE INVENTION 



Since the adoption of the synchronous optical network specification (SONET), a standard has been set 
for high-speed digital telecommunications (see American National Standards Institute, Inc. "Digital Hierarchy 
Optical Interface Rates and Formats Specification " standard TI.105 - 1988). Typically, such digital tele- 
20 communications combine a pseudo-random serial scrambling signal with the data stream so as to minimize 
the possibility of loss of clock signal which might otherwise result if the data stream comprised a large 
number of adjacent O's or Vs. However, due to the fact that the serial data stream may operate at 155 
megabits per second or higher, the serial PRG has to be implemented using high speed fabrication 
techniques, such as discrete emitter coupled logic (ECL) circuitry, ECL application specific integrated 
circuitry (ECL ASIC) or gallium arsenide (GaAs) circuitry, rather than the preferable CMOS circuitry which is 
less expensive to fabricate and operates at lower power than corresponding ECL or gallium arsenide 
circuitry. The additional fabrication costs and power requirements of ECL and GaAs circuitry also require 
more printed circuit board area in order to dissipate the additional heat, again making CMOS circuitry and 
especially CMOS application specific integrated circuitry (CMOS ASIC) preferable. 

Due to the fact that CMOS circuitry cannot typically operate at clock speeds greater than 50 megahertz, 
it is necessary that a technique be used to effectively reduce the clock frequency of the serial pseudo- 
random generator. The present invention describes such a technique and circuit which is operable for any 
serial PRG generating polynomial, as well as for any size parallel output word larger than the length of the 
equivalent serial shift register, representing the successive outputs from the serial PRG. 

In this manner, relatively low cost, low power consumption CMOS circuitry can be used to fabricate a 
parallel PRG which emulates the output of a serial PRG. 



SUMMARY OF THE INVENTION 



A parallel pseudo-random generator is described which emulates a serial pseudo-random generator 
which in turn operates upon a feedback arrangement wherein the next input value of the serial PRG is equal 
to the Exclusive-OR (XOR) combination of previous outputs of the serial PRG. For instance, in telecommuni- 
cations, a typical scrambling polynomial is 1 + x 6 + x 7 . This polynomial means that the next input value of 
the serial PRG is equal to the output of the sixth preceding value of the generator, exclusively ORed with 
the seventh preceding value of the generator. The output of the seventh preceding value of the generator is 
also typically exclusive ORed with the data to be scrambled. 

If the serial PRG has a clock rate of f SJ then the parallel PRG has a clock rate <f p ) of f s /W, where W is 
the number of outputs of the parallel PRG. 

The parallel PRG can be extended to any number of outputs (any size W) by choosing feedback paths 
which effectively emulate the serial PRG. The feedback paths are based upon the serial generating 
polynomial and the output size of the parallel PRG implementation. In two preferred embodiments of the 
present invention where W equals 8 and 16 respectively, corresponding numbers of D type flip-flops (FF) 
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are used with Exclusive OR (XOR) gates which provide the necessary feedback for determining the values 
of the next W outputs corresponding to the next W successive values of the simulated serial PRG. These 
two implementations are optimized using an optimization criterion set to the minimum number of exclusive 
OR gates for simulating the serial pseudo-random generator. 

OBJECTS OF THE INVENTION 



It is therefore a principal object of the present invention to provide a parallel pseudo-random scrambler 
circuit and method for simulating the output of a serial pseudo-random generator, such that W parallel 
outputs emulate W successive output values of the serial pseudo-random generator. 

Another object of the present invention is to provide a parallel PRG of the above description, wherein 
the value of W can be made arbitrarily large so that the resulting parallel clock frequency can be set 
arbitrarily low and therefore provide for implementation of a parallel PRG using CMOS fabrication 
techniques. 

A still further object of the present invention is to provide a parallel PRG of the above description 
incorporating D type flip-flops in association with exclusive OR gates for providing the necessary feedback 
from the W outputs so as to determine the next W outputs. 

A still further object of the present invention is to provide a parallel PRG of the above description which 
is implementable for any serial PRG generating polynomial. 

Other objects of the present invention will in part be obvious and will in part appear hereinafter. 

THE DRAWINGS 



For a fuller understanding of the nature and objects of the present invention, reference should be made 
to the following detailed description taken in conjunction with the following drawings, in which: 

Figure 1 is a diagram showing a serial pseudo-random generator incorporating use of D-type flip- 
flops connected as a P stage shift register, such that the next value generated is defined by the polynomial 
1 + X M + X* 

Figure 2 is a block diagram showing an 8 bit parallel PRG for emulating the serial PRG shown in 
Figure 5B. 

Figure 3 is a schematic diagram of the 8 bit parallel PRG shown in Figure 2. including clocking 
signals. 

Figure 4 is a block diagram of a 16 bit parallel PRG implementation of the serial PRG shown in 
Figure 1. 

Figure 5A is a diagrammatic representation for a W output parallel PRG implementation of a serial 
pseudo-random generator shown in Figure 5B, showing the feedback relationship between state m and 
state m-1. 

Figure 5B is a diagrammatic representation of a serial PRG similar to that shown in Figure 1 , wherein 
stage P and stage P-1 are the feedback values used to determine the next value of stage 1. 

Figure 6 is a transition matrix for the general solution of a parallel implementation of a serial pseudo- 
random generator corresponding to the parallel PRG shown in Figure 5A. 

Figure 7 is a diagram indicating the relationship between output (n) and the values of outputs (n + 6) 
and output (n+7) for the polynomial 1 + X 6 + X 7 . 

Figure 8 is a diagram indicating the relationship between output (n) and the values of outputs (n + 2), 
output (n + 5) and output (n + 9) for the polynomial 1 + X 2 + X 5 + X 9 . 

BEST MODE FOR CARRYING OUT THE INVENTION 



There has traditionally been a need to use serial pseudo-random generators for scrambling tele- 
communication information. As shown in Figure 1, a typical serial pseudo-random generator 20 (serial PRG) 
incorporates a plurality of stages arranged as a shift register 22 such that the value in each stage is 
transferred to the next stage until the last stage is encountered. The value in the last stage is typically 
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Exclusive ORed (XOR) with one bit of the telecommunication data stream with the result of the XOR 
operation actually transmitted in the telecommunication application. An Exclusive OR operation is defined 
such that if both inputs are logic 1 or logic 0, then the output is logic 0 and if the inputs are respectively 
logic 1 and logic 0, or vice-versa, then the output is logic 1 . A truth table representing an Exclusive OR 
5 operation is shown in Table 1. 

TAPLE 1 



Truth table for a two input EXCLUSIVE OR gate 
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addition with no carries) 



25 

The purpose of the pseudo-random generator in most telecommunication applications is to insure that 
regardless of the telecommunication bit stream pattern, the actual information transmitted will comprise 
approximately the same number of Vs and 0 T s. This result facilitates maintaining clocking information in the 
telecommunication bit stream which otherwise would be more difficult if, for instance, the telecommunication 
30 bit stream contained long, consecutive patterns of 1's or 0*s. Such scrambling is also useful for data 
encryption. 

Again referring to Figure 1 , it is seen that the operation of the serial pseudo-random generator can be 
defined by a polynomial typically of the nature 

1 + X M + ... + X p . This is known as the Characteristic Polynomial where w + n means the Exclusive OR 
35 operation (used in this manner throughout this document) 

The feedback equation associated with this characteristic polynomial is derived as follows: 

X° + X M + ...+ X p = 0 

Then, 

X c + X° + X M + ... + X p = 0 + X° 
to Since in general X + X = 0 and 0 + X = X where + is an Exclusive OR operation, 
then, 

X M +...+ X p = X° 

This last equation means that the next input value of the shift register is X M + ... +X P . 

For instance in the synchronous optical network (SONET) standard (also known as the American 

45 National Standards Institute (ANSI) standard T1 .105-1988), the polynomial is 1 + X s + X 7 . As seen in 
Figure 1 , this polynomial means that the value in shift register 6 is Exclusive ORed with the value in register 
7 with the result being the next value in stage #1 of the P stage shift register. Table 4 shows the values in 
the seven stages where the starting value for each of the seven shift registers is logic 1. This starting value 
is typically called a "seed". For the SONET standard, the seed is typically all Vs for a serial PRG. As is 

50 seen, the values generated for stage 1 successively move down the stages of the shift register. As noted 
above, the output from stage 7 is also used for Exclusive ORing with the serial telecommunication bit 
stream. 

The reason such a generator is called a pseudo-random generator is that the bit stream generated is 
always the same for the same starting seed and same polynomial. 
55 Although the SONET polynomial used Exclusive OR's stage 6 and stage 7, other polynomials may of 
course be used in which different stages of the serial shift registers are Exclusive ORed together. In fact, 
more than two stages can be Exclusive ORed if desired. 

Usually maximal length polynomials are used, that is, polynomials that repeat themselves after a 

4 
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maximum number of counts (clock cycles). For a maximal length polynomial the maximum number of 
counts is 2 n -1 for an nth order polynomial. For example, for a polynomial of degree equal to three, a 
maximal polynomial is 1 + X 2 + X 3 , while a non-maximal polynomial is 1 + X 1 + X 2 + X 3 . As seen in 
Tables 2 and 3, the maximal length polynomial repeats after seven outputs, while the non-maximal length 
s polynomial repeats after four outputs. 

The present invention is applicable witfi any serial polynomial, whether maximal or not. 
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Such a serial pseudo-random generator presents problems in integrated circuit implementation when 

40 the transmission rate of the telecommunication bit stream exceeds approximately 50 megabits per second. 
At speeds in excess of 50 megahertz per second, the fabrication of complimentary metal oxide silicon 
(CMOS) integrated circuitry becomes impractical. In fact CMOS fabrication at usable speeds exceeding 
approximately 75 megahertz is virtually impossible. As a result, for high transmission speeds such as those 
used in the SONET standard (such as 155 megabits per second), it is necessary if such a serial pseudo- 

45 random generator is to be used, that it be fabricated using emitter coupled logic (ECL) or gallium arsenide 
(GaAs) technology. Both these technologies have significant drawbacks as compared to CMOS technology 
in that they are typically more difficult to fabricate, and generate much more heat thereby requiring more 
printed circuit board area for placement of the integrated circuit components in order to dissipate the 
resulting heat, and cost more per logic gate. 

50 The present invention provides a general solution to the generation of high-speed pseudo-random bit 
patterns by providing a parallel pseudo-random generator having a plurality of parallel outputs whose values 
represent successive outputs of the serial pseudo-random generator. Such a parallel pseudo-random 
generator 24 may have any desired number of parallel outputs with the example shown in Figure 2 having 8 
outputs and that shown in Figure 4 having 16 outputs. The size of the parallel pseudo-random generator can 

55 be set to whatever value is best suited for a particular application, as long as the parallel word size is equal 
to or greater than the order of the scrambling polynomial. When using digital integrated circuitry the number 
of outputs generally has a value equal to a multiple of 2, such as 8 outputs, 16 outputs, etc. 

In the example shown in Figure 7, the pseudo-random generator comprises eight latches 26, which may 

5 
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be D type flip-flops, whose outputs (Q0 through Q7) represent eight successive output values of the 
emulated serial pseudo-random generator. Referring to Table 4 where the output of the serial pseudo- 
random generator is serial stage #7, it is seen that this seventh stage has logic value 1 for the first seven 
serial clock cycles (serial clock cycles 0-6) and has logic value 0 for the next clock cycle (serial clock 
5 cycle 7). Outputs Q7 through Q0 of the 8 bit parallel PRG therefore can represent these eight successive 
output values of stage 7 in the serial PRG as shown in Table 5. Thus it is seen in Table 5 that the Q0 
output represents the eighth sequential output of this serial PRG output stage 7, Q1 represents the seventh 
serial output of stage 7, and in similar fashion, down to Q7, 
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Table 4 Continued 

Serial Pseudo-Random Generator 

Corresponding to 1 + X 6 + x 7 generating polynomial 

Time Serial Stage # 

(equivalent aerial) 

(clock cycles) 12 3 4 5 6 7 
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50 1001010 

51 1100101 

52 1 1 1* 0 0 1 0 

53 1111001 

54 1111100 

55 0111110 
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Table S 
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55 which represents the first sequential output of stage 7 of this serial PRG. This pattern repeats for each new 
parallel output. 

As will be discussed in detail below, the next eight outputs of the parallel PRG from Q7 through QO 
have values 00000100. These values for Q7 through QO represent the next eight time sequential outputs of 

7 
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serial stage 7 as seen by comparing time outputs eight through fifteen of stage 7 presented in Table 4 with 
the parallel outputs for parallel clock cycle 1 (see Table 5). It is therefore seen that the initial (Zeroth) frame 
of the parallel PRG corresponds to the first eight sequential outputs (serial clock cycles) of stage 7 of the 
serial PRG, that frame 1 of the parallel PRG corresponds to the next eight sequential outputs of stage 7 
5 (serial clock cycles 8 through 15), etc. The initial parallel frame is the parallel seed input to the generator in 
order to start its operation in emulating the serial PRG which itself has a particular starting sequence or 
seed. 



w Analysis of the Parallel Implementation for Simulating a Serial PRG 

As seen in Figure 2, in addition to latches 26, the parallel PRG further incorporates a plurality of 
Exclusive OR gates 28 which combine various outputs of the latches for presentation as inputs to the 
latches for generating the next outputs on the latches. Figure 3 is a schematic diagram corresponding to 
15 Figure 2 showing additional logic circuitry for enabling the generator (AND gates 34), for loading the parallel 
seed (OR gates 36). and for presentation of a parallel clock signal 38. 

As seen in Figure 2, the inputs DO through D7 of the eight flip-flops, are presented with the values 
associated with functions FO through F7. These functions are defined by the equations presented in Table 
5A. 



TABLE 5A 
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It is further noted in Figure 2 that the outputs QO through Q7 of the parallel pseudo-random generator in 
turn are presented to a corresponding number of data stream output Exclusive OR gates 30 where the 

35 second input to each Exclusive OR gate is one bit of the serial data stream, such that the input to the Q7 
Exclusive OR output gate 30' is Exclusive ORed with the first bit of the serial data, output Q6 is Exclusive 
ORed with the next serial bit of data, etc., through QO which is Exclusive ORed with the eight bit of serial 
data. The output signals on output lines 32 therefore represent the scrambled output data which can then 
be converted back to a serial bit stream through use of an 8 bit multiplexer (not shown). 

4d it is readily seen in Figure 2, that if the parallel pseudo-random generator has a width of 8 (W = 8), that 
the frequency of the parallel operation is one-eighth that of the incoming serial bit stream since each 
parallel computation computes the next 8 outputs of this simulated serial pseudo-random generator as 
presented at outputs Q7 through QO respectively. 

45 

Determination of the Parallel Output Exclusive OR combinations 

As will be presented more fully below, the determination of the Exclusive OR gate arrangement for 
presentation as an input to each of the parallel pseudo-random generator latch is determined in a manner 

50 so as to emulate the serial pseudo-random generator output bit stream. Although a particular Exclusive OR 
gate arrangement is shown in Figure 2, there are in fact many implementations which are possible. The 
present invention is particularly advantageous when the minimum number of Exclusive OR gates are used 
for each input. This arrangement minimizes the requirements for serial gates and consequently minimizes 
the gate delays associated with each serial gate. 

55 It has been experimentally found and mathematically verified as presented hereinafter in a mathematical 
analysis by inventor G. Roger entitled "Parallel Pseudo-Random Generator, Mathematical Analysis", that for 
any serial pseudo-random generator polynomials, there exists a solution by which Exclusive OR gates can 
be used to implement a parallel pseudo-random generator provided that the width of the parallel PRG is at 

8 
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least equal to the maximum shift register stage used to define the serial PRG polynomial. 

For the polynomial presented above with regard to Figure 1, that is, wherein the next input to stage #1 
is equal to the Exclusive OR output of stage 6 and 7, it is seen that this relationship can be defined 
generally as follows: 

5 (1) Q(n)s Q(n + 6) + Q(n + 7); where "n" is any stage of the serial PRG. Figure 7 shows a graphical 
representation of this relationship. 

Again, referring to Table 4, it is seen that stages 6 and 7 for clock cycle 0 both have a logic 1 value. 
Consequently, the next value for stage 1 is equal to 0 (1 +1 =0, see Table 1). This result is analogous to 
the above formula where n equals 0 (Q(0) becomes Q(1) after the next serial clock cycle, and in general Q- 
70 (n-1) becomes Q(n) after the next serial clock cycle). 

In order to determine the next 8 bits of the emulating eight output parallel pseudo-random generator, it 
is observed that the next generated bit of the serial PRG will become, after eight serial clock cycles, the 
next value for output Q7 of the parallel PRG (see Table 6 where Q-1 become Q0 after one serial clock 
cycle; which becomes Q7 after seven additional serial clock cycles; where these eight serial clock cycles 
75 are equivalent to one parallel word clock cycle). Thus for the eight bit parallel PRG implementation, the next 
value for Q7 is equal to Q-1 which is equal to the Exclusive OR of Q5 and Q6, that is: 
Next Q7 = F7 = Q5 + Q6. 

Using this same rationale, it is seen that the next value of Q6 through Q2 can be defined as follows: 



Next Q6 




F6 




Q4 


+ 


Q5 


Next Q5 




F5 




Q3 


+ 


Q4 


Next Q4 




F4 




Q2 


+ 


Q3 


Next Q3 




F3 




Q1 


+ 


Q2 


Next 02 




F2 




Q0 


+ Q1 



The evaluation of Q1 can also best be understood by reference to Table 6 and Figure 7. 

Table $ 
parallel output values for 
two eight bit words 
(from n - -8 to n - 7) 



Q-8 Q-7 Q-6 Q-5 Q-4 Q-3 Q-2 Q-1 | 4 QO Ql Q2 Q3 Q4 Q5 Q6 Q7 
next 8 bit word present 8 bit word 



35 Thus the next value for Q1 (which equals F1) is equal to the value of Q-7. 
Next Q1 = F1 = Q-7 
from equation (1): 

F1 = Q(-7 + 6) + Q(-7 + 7) « Q-1 + QO (n = -7) 

but Q-1 = Q5 + Q6 (n = -1) (using equation 1 again) 
40 therefore: 

Next Q1 = F1 = Q-1 + QO = Q5 + Q6 + QO 

However, it is also seen from equation (1 ) that the present value of QO is equal to the present value of 

Q6 + Q7 (QO = Q6 + Q7), and thus 

(2) Next Q1 = F1 = Q5 + Q6 + Q6 + Q7 
45 Since the Exclusive OR of any logic value with itself is equal to 0 (see Table 1 above), equation (2) can 

be rewritten as follows: 

Next Q1 = F1 = Q5 + Q7 

Using the same rationale, it is readily seen that the next value of QO is defined as follows: 

Next QO » F0 = Q4 + Q6 

50 Therefore the essence of the procedure for determining the Exclusive OR gate arrangement is to . 
determine through the serial generating polynomial, the inter-relationship between the serial stages. Since 
the parallel relationship merely displays a plurality of serial stages at the same time, then the serial 
polynomial is used to compute the next parallel output for each of the parallel outputs after W serial clock 
cycles, where W is equal to the width (i.e. number) of parallel outputs. Since only the present values of the 

55 parallel output stages are available for computing the next values of these same stages, if an output value is 
required from one or more of the next outputs (next word as shown in Table 6) of the parallel PRG, then the 
serial polynomial is again used for that particular output to determine the present outputs which represent 
that next output value. This procedure can be used with any parallel word size and for any serial generating 

9 



EP 0 397 079 A2 



polynomial. 

Referring to Table 7, it is seen that for a 16 bit parallel PRG implementation, the next value of Q15 
simply equal to Q-1 (that is the serial output 16 serial clock pulses later) and thus- 
Next Q15 = F15 = Q-1 = Q5 + Q6. (see equation (1) for n = -1) 
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This analysis holds for the next values of Q14 through Q10. that is 
Next Q14 = F14 = Q4 + Q5 
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Next Q13 = F13 = Q3 + Q4 

Next Q12 = F12 = Q2 + Q3 

Next Q11 = F11 = Q1 + Q2 

Next Q10 = F10 = Q0 + Q1 
s it is seen that the next value of Q9 equals F9 which equals Q-1 + QO. However. Q -1 Is simply equal to 

Q5 + Q6 and thus 

Next Q9 = F9 = Q5 + Q6 + QO. 

The present value of QO is, by equation (1), equal to Q6 + Q7 and thus: 

Next Q9 = F9 = Q5 + Q6 + Q6 + Q7. 
10 Next Q9 * Q5 + Q7 

Similarly, for the next values of Q8 through Q4 are: Next Q8 = F8 = Q4 + Q6 

Next Q7 = F7 = Q3 + Q5 

Next Q6 = F6 = Q2 + Q4 

Next Q5 = F5 = Q1 + Q3 
15 Next Q4 = F4 = QO + Q2 

The next value of Q3 is equal to Q-13. Using equation (1) above, we have the following: 

Next Q3 = F3 = Q-13 

Q-13 = Q-7 + Q-6 (n = -13) 

Q-13 = (Q-1 + QO) + (QO + Q1) 
20 Q-13 = Q-1 + Q1 

Q-13 = (Q5 + Q6) + Q1 

(2a) Q-13 = (Q5 + Q6) + (Q7 + Q8) 

but also from equation (1 ): 

Q5 = Q11 + Q12 
25 Q6 = Q12 + Q13 

Q7 = Q13 + Q14 

Q8 = Q14 + Q15 

so therefore, 

Q-13 = (Q11 + Q12) + (Q12 + Q13) + (Q13 + Q14) + (Q14 + Q15) and therefore: 
30 (2b) Q3 = Q11 + Q15. 
Similarly, 

Next Q2 = F2 = Q10 + Q14, 
Next Q1 = F1 = Q9 + Q13, and 
Next QO = F0 = Q8 + Q12. 
35 The Exclusive OR implementation for the 16 bit parallel PRG is shown in Figure 4 corresponding to the 
above analysis. 

It is therefore seen that the next Q13 can be defined by a plurality of Exclusive OR operations such as 
shown by equations (2a) and (2b). In general such multiple representations can be shown for the outputs. 
One optimization criterion can be to use the minimum number of gate inputs, which is shown by equation 
40 (2b) for output Q3. 

The above analysis can be used for any width parallel PRG provided that the width of the parallel PRG 
is at least equal to the maximum number of serial stages used in the feedback arrangement for the 
emulated serial PRG. In the example above, where the serial polynomial uses stages 6 and 7 to compute 
the next input stage, the value of P equals 7 and thus, the width of the parallel PRG must at least equal 7, 
45 although it may be any size greater. 

Furthermore, although the serial polynomial was equal to the Exclusive OR of two serial stages, the 
present invention is applicable to any serial polynomial, regardless of the number of serial stages 
Exclusively ORed used to compute the next input 

An example of a more general serial pseudo-random generator is the following polynomials: 
so Next serial input = X 2 + X s + X 9 . 

That is. the characteristic polynomial is 1 + X2 + X5 + X9. 

This polynomial is non-maximal (see Table 2 and 3 above) and is presented to demonstrate that the 
parallel PRG implementing methodology is general in application. 

Figure 8 diagrammatically shows this serial pseudo-random generator in terms regarding stage n f such 

55 that 

(3) Q(n)sQ(n + 2) + Q(n + 5) + Q(n + 9). 

Table 8 shows the serial stage values for the nine stages comprising the serial pseudo-random 
generator corresponding to this polynomial over 36 clock cycles (clock cycles 0 - 35). 
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By using the relationship in equation (3) the values for the next QO through Q8 outputs of a parallel 
pseudo-random generator with width W = 9 emulating the serial pseudo-random generator shown in Figure 
8 are as follows: 

Next Q8 = F8 = CM = Q1 + Q4 + Q8 

Next Q7 = F7 = Q-2 = QO + Q3 + Q7 

Next Q6 = F6 = Q-3 = Q-1 + Q2 + Q6 
= Q1 + Q4 + Q8 + Q2 + Q6 
= Q1 + Q2 + Q4 + Q6 + Q8 

Next Q5 = F5 = Q-4 = Q-2 + Q1 + Q5 
= QO + Q3 + Q7 + Q1 + Q5 
= QO + Q1 + Q3 + Q5 + Q7 

Next Q4 = F4 = Q-5 = Q-3 + QO + Q4 
= Q-1 + Q2 + Q6 + QO + Q4 
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= Q1 + Q4 + Q8 + Q2 + Q6 + Q0 + Q4 
= Q1 + Q8 + Q2 + Q6 + QO 
= QO + Q1 + Q2 + Q6 + Q8 

Next Q3 = F3 = Q-6 = Q-4 + Q-1 + Q3 
= (Q-2 + Q1 + Q5) + (Q1 + Q4 + Q8) + Q3 
= (QO + 03 + Q7) + Q1 + Q5 + Q1 + Q4 + Q8 + Q3 
= QO + Q7 + Q5 + Q4 + Q8 
= QO + Q4 + Q5 + Q7 + Q8 

Next Q2 = F2 = Q-7 = Q-5 + Q-2 + Q2 
= (Q-3 + QO + Q4) + (QO + Q3 + Q7) + Q2 
= ((Q-1 + Q2 + Q6) + QO + Q4) + (QO + Q3 + Q7) + Q2 
= (((Q1 + Q4 + Q8) + Q2 + Q6) + QO + Q4) + (QO + Q3 + Q7) + Q2 
= Q1 + Q8 + Q6 + Q3 + Q7 
= Q1 + Q3 + Q6 + Q7 + Q8 

Next Q1 = F1 = Q-8 = Q-6 + Q-3 + Q1 
= (QO + Q7 + Q5 + Q4 + Q8) + (Q1 + Q4 + Q8 + Q2 + Q6) + Q1 
= QO + Q7 + Q5 + Q2 + Q6 
= QO + Q2 + Q5 + Q6 + Q7 

Next QO = FO = Q-9 = Q-7 + Q-4 + QO 
= (Q1 + Q8 + Q6 + Q3 + Q7) + (QO + Q3 + Q7 + Q1 + Q5) + QO 
= Q8 + Q6 + Q5 
= Q5 + Q6 + Q8 

Table 9 shows the output values of the parallel pseudo-random generator for four parallel clock cycles 
corresponding to serial clock cycles 0 through 35. It is seen that these outputs correspond to the serial 
pseudo-random generator output stage 9 for the first 36 serial clock cycles. 

Table 9 

Parallel Pseudo-Random Generator 
' (width - 9 bits) emulating 
serial PRG corresponding to 

1+ X 2 + X 5 + X 9 polynomial 
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It is observed from the foregoing that as long as the width of the pseudo-random generator is at least 
equal to the number of stages used in the serial pseudo-random generator, the parallel pseudo-random 
generator is implementable. It is further seen that the minimum number of Exclusive OR gates necessary 
for implementing the parallel PRG is not necessarily equal to the number of Exclusive OR gates used in the 
corresponding serial PRG, at least when the parallel PRG has a width equal to the serial PRG. 

The following mathematical analysis proves that there always exists a parallel PRG emulation of the 
serial PRG when the width of the parallel PRG is at least equal to the number of stages needed to 
implement the serial PRG polynomial. 



Parallel Pseudo-Random Operator 
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Mathematical Analysis 



1. Introduction 

10 A parallel pseudo-random generator Is analyzed to replace a classical serial PRG generator built with 
shift registers. Both the parallel and serial generators are represented by the schematics diagrams shown in 
Figure 5A and 5B respectively. 

In the classical solution, the signals issued from several stages of a p-stage shift register are added 
together by Exclusive OR (XOR) gates, and the input of the register is fed with the resulting signal, creating 

rs a feedback. The equation between the successive values of the signal is: 

(4) Sn = Al Sn-1 + A2 Sn-2 + + Ap Sn-p 

in which ' + ' is used for XOR. or modulo 2 addition, and AI,...AL..Ap are 1 if the output of the stage i is 
connected, or 0 if not. This is an equation whose coefficients are in the field of integers modulo 2 'FfO.D'. 
The *Z transform* of the signals leads to: 
20 (5) S(Z) = Al r S(Z) + A2 Z 2 S(Z) + ... + Ap Z P S(Z), or 

(6) P(Z). S(Z) = 0, with 

(7) P(Z) = Z° + Ai Z 1 + A 2 Z 2 + ... + ApZ p Equations 5 and 6 are equivalent 
P(Z) Is the characteristic polynomial of S and may be considered a 'generator' of S. 

If the polynomial P(Z) is 'irreducible and primitive' (is not a product of polynomials of smaller degrees 
25 with coefficients in F(0.1)), and has a primitive root of Z q + 1 = 0 (g = 2 P -1), the sequence generated by 
the system will be a pseudo-random generator of period 2 P -1. 

The parallel generator consists of a multi-output latch (e.g. a plurality of flip-flops), the input signals of 
which are computed by a network of XOR gates, this network being fed by the output signals of the latch. 

30 

2. Preliminary Remarks 

These remarks may be useful to the reader unfamiliar with the methods of digital signal processing. 

1) The field of 'integers modulo 2' contains only two elements, namely 0 and 1, with two operations; 
35 multiplication (AND) and addition (Exclusive OR, or XOR) such that: 

0x0 = 0,0X1 =1X0,1X1 = 1 and 
0 + 0 = 0, 0 + 1 =1 + 0 = 1,1+1=0 

Polynomials with their coefficients in this field have properties such that: 
P(Z) = Q(Z) is equivalent to P(Z) + Q(Z) = 0. or 
40 (1 + Z) 2 = 1 + Z 2 (because 2 = 0). 

2) We use a r Z transform* such that if S n = S(nT) (T is a time interval) 
S (Z) = E S n Z n Z is a 'lag operator', because 

ZS(Z)+ IS n Z n * 1 = ES n . t Z n , 
which is the Z transform of S(t-T). 
45 3) If we look for the solutions of: 

(5) S n = At Sn., + A 2 S n - 2 + + A p S n .p , it is usual to let: 

S n = C a" 1 *, C being a constant. 

The equation becomes: 

a" n = Av a*"* 1 + A 2 a n " 2 + ...+ A p a n * p , or 

so a 0 + An a + A 2 a 2 + + A p a p = 0 

and we see that "a" must be a root of P(Z) = 0. There are p roots of this equation. 

These roots generally cannot be expressed with 0 or 1 . but a general solution of (5) will be a linear 
combination of the successive powers of them: 

S n = a; n + a 2 * + .... + a D n , which is a symmetrical function of the roots of P(Z) t therefore a function of the 
55 coefficients of P(Z), which are equal to 0 or 1 , if these coefficients are on the field of the integers modulo2. 

4) It is easy to see that the shift register of the serial PNS generator which contains p bits may 
contain at most 2 P words, including the word 0, 0, 0... which generates a null sequence. Therefore, there 
may exist at most 2 P - 1 non null words and the period of the sequence is at most 2 P - 1 . That period is 
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obtained with particular polynomials called 'irreducible and primitive 1 polynomials. 
5) If w a w is a root of P(Z) = 0. 

P(a) = a 0 + Ai a 1 + A 2 a 2 + + A p a p = 

P 2 (a) = a 0 + Ai a 2 + Aza 4 + + A p a 2p , 

5 

because A^ -» A^ 
p-1 

70 

and a 2 is another root of P(Z) = 0. as are a...a 2 , which are the p roots of the equation. The next one, with 
exponent 2 P is equal to "a" because 

P 



20 

3. THE PARALLEL GENERATOR. 

The 'Z equations 1 of the parallel generator of Figure 5 A are: 



25 



-0 t _ -N -N+1 A Z 2N-1 

Z A 00 2 + A 01 Z + A O f N-l Z 



30 



Z l - A i0 Z* * A A1 Z N+1 ...... A i>M _ 1 Z 2 "' 1 



40 



K-l N K+l 2N-1 

2 - A N-1,0 2 + *N-1,1 2 + •'♦• + A N-1,N-1 2 



The matrix (the elements of which are A^j) is the transition matrix between two successive states m-1 
45 and m of the latch. These coefficients are 1 or 0, following output j is linked or not to input.!, generally 
through XOR circuits. For example, equation i corresponds to: 



50 



S n-i- < A i0 W XO * < A iI S n-H-l> * OR < A i2 S n-N-2>-** 
Equation i may be written: 

z 1 - r^ 1 h z N+ 3 or z° - z"" 1 x*' 1 a z* -z^ 1 R (Z) 
j-o ij j-o ij i 
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Ri is a polynomial, the coefficients of which are the elements of row r of the transition matrix. 
Tj (2) = Z° + Z Ni Rj must be such that T^Z) .S(Z) =0. 

We know that P(Z) .S(Z) = 0 f therefore, if T(Z) is a multiple of P, for example T(Z) = P(Z) .G(Z), then 
T(Z) .S(Z) = P(Z) .Q(Z) .S(Z) = Q(Z) .[P(Z) .Z(Z}] = 0 
s (This result may be obtained by considering that successive values of S are combinations of powers of the 
roots of P(Z) = 0, which implies that T(Z) = 0 must have at least the same roots as P(Z) = 0). 

Suppose now that we consider a polynomial P = AO + AiZ + A 2 Z 2 + A 3 Z 3 . 
(we take a particular example, easy to understand, but the derivation is general). 
The sequence generated by P(Z) is such that: 
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Let Q(Z) = Bo + B*Z + B2Z 2 

then 

T(Z) = P(Z) .Q(Z) = A c Bo+ (AoB, + A,B 0 ) Z + (A 0 B 2 + AiB, + A 2 B 0 ) Z 2 + (AiB 2 + A 2 Bi + A 3 B 0 ) Z 3 
+ (A 2 B 2 + A3B1) Z* + A3B2Z 5 , 
25 orT(Z) = C0+C1Z +C2Z 2 + C 3 Z 3 + C4Z 4 + CsZ 5 
If we compute, S being a sequence: 
CoS'n.5 + CiS' n . 4 + ... +CsS'„ = A 0 B 0 S'„* 5 + (Ao B1 +A1B0) S'„ + 4 + ..., 
we get: 

B 0 ( Ao S n/ K 5 + Ai S'„* 4 + A 2 S„«. 3 + A 3 S' n *2 ) + + Bi (Ao S'„* 4 + Ai S' n + 3 + A 2 s'„* 2 +A 3 S' n *i ) + + 
30 B 2 ( Ao S'„ +3 + A1 S'n :2 + A 2 S'„*i +A 3 S'„ ) 

which is equal to 0 if s' is the sequence generated by P(Z). 

Furthermore, we see that if S„, S n *t, S' n . 3 , s' n *4 are successive values of S and if the computed sum 
is null. S n-5 »s the following sample of S. 

We conclude that if T(Z) is a multiple of P(Z), 
35 1) T(Z) S(Z) = 0, which was foreseen, and 

2) T(Z) generates the sequence S, if it is fed with a 'good seed 1 . This result means a series of 
samples of the sequence S (with another seed, it could generate a sequence generated by Q). 

T(Z) is such that only its first coefficient (always equal to 1), and its N last coefficients (those of Rt). 
may be different than zero. Therefore, the needed seed is limited to the samples contained in the latch. At 
40 the starting time, it is necessary that the latch be loaded with a section of the sequence generated by P(Z). 
Every polynomial T t (Z) generates a sample of the following series of bits of the latch, and for each clock 
time, the series of N bits located in the left part of the figure (state m) becomes the" series located in the 
right part (state m-1 ), etc. 

Figure 6 gives a picture of the coefficients of polynomials Tj (Z). The, coefficients equal to 1 are noted 
45 'X\ the others being equal to 0. The coefficients of R lf which are the elements of the transition matrix, are 
inside a parallelogram, and we see that our problem is to find polynomials which: 

1) are multiples of P(Z), or generators of the sequence generated by P(Z) (which is equivalent), 

2) have their coefficients, other that the first one, included in the parallelogram, 

3) have a minimum of terms equal to 1, in order to yield the simplest implementation. 

50 The 'Bezout's Relation* (see Table 10) allows one to compute polynomials which fulfill the first two 
conditions, but not always the third condition. 

A very simple method to find the 'good* polynomials is to 'try and see 1 : for each line of the matrix, 
polynomials with two non-null coefficients included in the parallelogram. These polynomials are tested as 
generators of the sequence S and the first one being found is used and, if possible, reused for the next 
55 rows of the matrix. If searching with two coefficients fails, we look for three coefficients or more, etc. 

Even though the polynomial is simple in nature, it may required tens of seconds of computing time (20 
seconds with a D.E.C. VAX 8600 Computer for a polynomial of the 12th degree and N = 32, but virtually 
immediate results for a polynomial of the 7th degree and N = 8 or 16) Note N here is equivalent to W; i.e., 
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number of parallel outputs. 



4. Alternate Derivation 



a) The Characteristic Equations of the Parallel System 

Such a parallel generator may be represented by two successive states of the latch, linked by a 
transition matrix, the elements of which are 0 or 1, 1 meaning an XOR operation. So. each signal of the 
second state depends on, the signals of the first one by: 

(8) 7} = L j Bij Z' + ft with 0 £ i £ N-1 and 0 £ j < N-1 , 
Bij = 0 or 1 . 

N is the number of bits contained in the latch and a series of N (N is equivalent to W as presented 
earlier) bits of the PRG are delivered for each parallel clock cycle instead of one bit with the shift register 
solution per serial clock cycle. Letting k = N-i, we may replace (8) by: 

(9) Z Nk = Z N N I j Bkj 7) = Z N Rk t where 

Rk = E k Bkj Z j is a polynomial of the (N-1)th degree with Bkj coefficients and the kth (or N-i)th) row of the 
matrix. 

We replace equation (8) by: 

(10) Z° = 1 = Z k Rk 
or: 

(11) Sk(Z) = 1 + Z k Rk = 0 

The equation Sk(Z) = 0 is a characteristic equation which must be fulfilled for every row k of the 
transition matrix, with 1 < k < N. 

b) Properties of the Characteristic Equation 

The successive powers of the p roots of P(Z) are solutions of equation (4). The signals Sn are linear 
combinations of these powers, and are symmetrical functions of the roots of P(Z). 

Therefore, to generate the same signals as the serial shift registers, the roots of Sk(Z) must include the 
roots of P(Z), and: 

(12) SK (Z) = P(Z) Qk(Z) 

where Qk(Z) is a polynomial which must have at least a Z° term, since Sk and P have such a term. 
Does the converse hold? 

If equation (11) : Sk (Z) = P(Z) Gk(Z), then Sk(Z)=0 is true not only for the roots of P(Z). but also for 
the roots of Qk (Z). There could be a problem with what are known as parasitic roots, but by choosing a 
•good seed 1 (that is. a segment of the good PRG), we avoid introduction of such parasitic roots, as it may 
be proved, considering the successive values of Sn : if Sn. Sn-1, etc. are contained in the latch and are a 
segment of the PRG. Sn + 1 will be a bit of the same PRG, and if Sn + 1. Sn, Sn-1, ...are a segment of the 
PRG, then Sn +2 will be a bit of the PRG etc. 

So. every polynomial Sk(Z) must be a multiple of P(Z). starting with Z°. Such a polynomial has other 
terms only between Z k and Z k * NO (the terms of Rk), and conversely, such a polynomial is convenient for a 
parallel generator. There may be several equivalent expressions of Sk(Z) for the same value of k. 

There may exist different expressions of Sk. For example, the following two polynoms may be valid: 

(13) Sk1 (Z) = 1 +Z k Rk1 = P(Z) Qk1(Z) 

(14) Sk2(Z) = 1 + Z k Rk2 = P(Z) Qk2(Z) 

The only condition is that degrees of Rk1 and Rk2 both be less than N. 
By subtraction, we obtain: 

(15) Sk1 - Sk2 = Z k ( Rk1 - Rk2) ■ P(Z) ( Qk1 - Qk2) 

First, the polynomial (Sk1 - Sk2) is divided by P(Z). Sk1 and Sk2 are said to be 'congruent modulo P- 
(Z)\ It means that Sk2 may be obtained by replacing terms of Sk1 , with reference to P(Z). 

For example, if P(Z)= 1 + Z* + Z 7 , we may, in Sk1, replace Z by Z 7 + Z 8 . because if 1 + Z 6 + Z 7 
= 0. Z + Z 7 + Z 8 = 0 also. 

Second, as Z k cannot divide P(Z), which is a prima, without a null root, it divides (Qk1 -Qk2). So (RK1 - 
RK2) is also a multiple of P(Z) and Rkl and Rk2 are also congruent modulo P(Z). 

Thus there may exist, for the same row of the matrix, several polynomials Rk(Z), 'congruent modulo P 
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(Z) 1 , (with their degrees less than N), which give equivalent characteristic polynomials Sk(Z). 
Figure 6 shows several aspects of the problem: 

1) There are two coordinate systems to take into account: one for the polynomials Sk, and one for the 
polynomials Rk, which must lie inside or on the edges of a parallelogram (from k = 1 to k = 8) and the 

5 allowed positions for 1 coefficients of the Rk are marked 7\ 

2) In the example chosen, there 'are two interesting multiples of P(Z); namely, 1 + Z s + Z 7 and 1 + 
Z 12 + Z 14 , 

and we see that the non-constant terms of Sk may be non-constant terms of one of these multiples, if all of 
them are in the allowed domain. 
10 3) Several polynomials Sk may be identical (S1 to S6 for example), and the corresponding Rk will 

differ only from a translation of terms. 



c) The Bezout's Relation 

15 

(16) Sk(2) = 1 + Z k RK = P(Z) Qk(Z), or 

(17) 1 = Z k Rk + P(Z) Qk(Z) 

We recognize the BEZOUTs relation (see 'Relation of BEZOUT later, Table 10): 

if P(Z) and Q(Z) are two polynomials, their greatest common divider (GCD, or Hcf- highest common 
20 factor -) may be expressed as: 

(18) HCF (P,Q) = A(Z) P(Z) + B(Z) Q(Z). 

A and B can be found with a very simple algorithm, (derived from the Euclidean algorithm) and the 
degree of B is smaller than the degree of P. 

Z k and P(Z) have no common factor (P is irreducible) so their HCF is 1, and equation (17) is the 
25 BEZOUTs relation. 

So, for every value of k. we are able to determine a polynomial Rk, the degree of which is smaller than 
p, degree of P(Z). Taking k = 1 in (17), we see that since 1 = Z R + P Q, and since the degree of P Q is at 
least equal to p, and degree of R is at most equal to p, the unique possibility is degree of R = p-1, with 
degree of Q being 0. So, at least one of the polynomials Rk has p terms and the transition matrix must have 
30 at least p columns. 

Therefore: 

1) N must be at least equal to p (see discussion above concerning observed relationship between W 
(the same as N here) and P). 

2) For N greater or equal to P, there is at least one solution to our problem. 

35 That solution, in general, will not be optimum because we are typically looking for a minimum of XOR 
circuits. But if N > p-1, we shall have a way of improving the solution by searching polynomials Rk 
congruent modulo P(Z) with those given by the BEZOUTs relation, with the further condition that their 
degrees be less than N. 

40 

d) The 'Heuristic Solutions'. 

The heuristic solutions consist of searching systematically for the multiples of P(Z) having two. three, or 
more non-constant terms. If a two coefficient solution is found for the row k, it is used as much as possible 
45 for k + 1,... If three or more coefficients are needed, they are used only for the row k, because we may hope 
that the following row will accept less coefficients, and we start again with two coefficients. To test we divide 
by P(Z) and look to see if the remainder is the null polynomial. It may become expansive in computing time 
for high values of P and N, but it leads sooner or later to an optimum solution (there may be several 
solutions). We take the first one we find, at least at the present time. 
5a Another way to test an Sk polynomial is to verify directly that the polynomial is able to generate the 
pseudo-random sequence generated by the given characteristic polynomial. This method is used in the 
program 'GSPA-E' (see Table 11 and the TEST portion of subroutine POLYANCOEFX). 

Of course, other strategies are possible, depending on the objective. For example, it could be better to 
compute a table of the polynomials Sk able to fulfill the conditions, and pick among them to build the 
55 matrix. i 

Of course, such multiples as (1 + Z 5 + Z 7 ) 2 are evident as good. 



18 



EP 0 397 079 A2 



In Summary 

The problem is to know the multiples of P having a minimum of coefficients and pick among them 
those whose non-constant terms fall in the range of the Rk polynomials. 



e) The Program 

There are four parts to the program: 
ro 1) Initialization and input of the data: 

- degree of the polynomial P(Z), p 

- coefficients of P other than AO and Ap (which are always 1) 

- number of bits of the latch, N 

The PRG sequence corresponding to P(Z) is generated (Seq1), in order to: 
75 - make sure that P is a 'good 1 polynomial 

- prepare a good 'seed 1 for the test of the parallel system 

- have a reference to test the parallel system. 

2) Computation of the matrix elements 

3) Publication of results 

20 - a table of the coefficients of the matrix 

- a drawing of the matrix, if N is no more than 32. 

4) Verification 

- a sequence (Seq2) is generated and compared to Seq1. 

A file of subroutines is used with these programs. It contains all of the operations in modulo-2 algebra 
25 needed for our purpose. 

Table 12 is a sample terminal listing from execution of the GPSA-E program. 

Table 13 contains several printouts for 8, 16. 24, 32, and 64 bit parallel word widths of the SONET 
polynomial (1 + X s + X 7 ). 

30 
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Table, 10 

•Relation of Bezout' 

We want to obtain the Highest Common Factor (HCF) of two integers, 
a and b, or of two polynomials. The algorithm is similar for both. 

First , we divide a by b: 
a - QO b + Rl 0 i Rl< b 

The HCF, which divides a and b, divides Rl. We divide now b by Rl, 
etc . * . , 
b - Ql Rl + R2 Q£ R2 <R1 

Rl - Q2 R2 + R3 QS R3 <R2 

. R is smaller and smaller. Therefore, 

it happens: 

Rn-2 « Qn-1 + Rn-1 + Rn 

(Rn is the HCF) 

Rn-1 - Qn Rn + Rn+1 with Rn+1-0 

and Rn, which divides Rn-1, divides Rn-2 , a, and b. 

Rn is the HCF of a and b. This is the Euclidean Algorithm. 

Now consider the sequence of the successive remainders: 
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5 



10 



15 



20 



25 



35 



Rl- a -QO b - Al a + Bl b, with Al- 1 and Bl- -Q0 

R2- b -Ql Rl - A2 a -f B2 b, A2~ QX Al, B2- -Ql Bl 

R3- Rl -Q2 R2 - A3 a + B3 b, A3- Al - Q2 A2 , B3- Bl -Q2 B2 



Rn- An a + Bn b ,An - An-2 -Qn-1 An-1, Bn- Bn-2 - Qn-1 Bn-1 

Therefore, An and Bn are obtained from Al at Bl. If ve sat: 

A- 1-1, B-1-0, (-1 is a subscript), and AO - 0 and BO - 1, tha algorithm 

which yields Rn, An and Bn, starting fro» subscript 1, is simple to 
implement. m 

If Rn is the HCF, we obtain the 'relation of Bezout 1 ; 

IICF (a,b) - A a + A b 



If a and b are coprime ('premiers entre eux 1 in f rench) , 
Rn - HCF(a,b) - 1. If a and b are polynomials, Rn is a constant, 
30 which, if the coefficients are in F(0,1), is equal to l, and we obtain: 

! 1 

I 1- a(Z) A(Z) + b(Z) B(Z) 1 
j * -I 



40 



It may be seen, or computed, that the degree of An is the degree of the 

product; Ql .Q2 • , , .Qn-1 , and for Dn, it is the degree of Q0.Q1 Qn-1. 

Furthermore: ('A means degree of polynomial A) 
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10 



75 



20 



25 



30 



35 



40 



45 



50 



•Q0« - a - # b 

•Q2- *R1 - -R2 

♦Qn-1 - «Rn-2 - •Rn-l 
•Qn - *Rn-l - *Rn 

Therefore, • (qo . qi . Q2...Qn-l) - •* - -Rn-l 

If we suppose that a and b are coprime, and that Rn is their HCF, 
•Rn - o and -Rn-l ie at least 1. therefore, «Bn is smaller than -a, 
and by similar reasoning, "An is smaller than *b. 

End of Table 10 



55 
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TABLE 11 



C k * PROGRAM GPSA £ 

C * COMPUTAT lOHTO^'lIE TRANSITION MATRIX or A PARALLEL GENERATOR • 

c or PSEwOORANDOM SEC/UENCE • 

c ...*.(C. ROGER. LABORATOIRE5 DE MARCOUSS I S . OCTOB ER 9 1988)« 

C • rORTRAN-7 7 VAX VMS 

<;*•• DECLARATIONS 

c 

CHARACTER* 3 NON,*EP 
CHARACTERS NOM4 
C1IARACTER*75 LIGNE 
INTEGER DEC 

PARAMETER ( DEG- 1 30 . NB RN-6 4 , NBB ITS- 1 0000 ) INBBITS POUR DEGRE<-12 
C •* POLYNOMIALS LIMITED TO DEGREE 1 2 8 . PARAMETERS MAY BE AJUSTED AS NEEDED 
C NBBITS IS AT LEAST TWICE THE NUMOER Of BITS Or THE P.N.S. 

C P IS THE DEGREE Or THE CHARACTERISTIC POLYNOMIAL GIVEN TOR THE SERIES 

C GENERATOR . P 2 IS THE CHARACTERISTIC POLYNOMIAL ITSELF . 

I HTEGER P . P i i DEG ) , P 2 ( DEC ) , 0 ( D EG ) , R « DEC ) 

INTEGER RN ( NBRN , DEG ) I THE TRANSITION MATRIX 

INTEGER A(DEG) 

INTEGER B(DEG) 

INTEGER SE01 (NBBITS) , S EQ2 ( NBB ITS ) 
C SEOI AND SE02 ARE TWO P.N.S. 
C 

C INITIALISATIONS AND INPUT OT PARAMETERS 

C *• •»*****»«* *» 

1 1 DEG- DEG 

IUNIT1-6 1 SECONDARY RESULTS (FOR VER t T ICATIONS ) 

IUNIT2-6 !THE SCREEN , 

IUNIT3-5 1THE KEABOARD 

IUNIT4-10 I IMPORTANT RESULTS (ON THE TILE GPSA. DAT) 

OPEN ( UN IT* 1 0 , F I LE- ' GPSA . DAT* , STATUS— ' HEW* ) 

WRITE (IUNXT2,J01) 

WRITE (IUNIT4,J01) 
301 fORMAT <//.' PARALLEL PSEUDONOISE SEQUENCIES GENERATOR. * ,/ $ 

C COMPUTATION Or THE TRANSITION MATRIX. ( PROGRAM GPSA)',/, 

C ' GEORGES ROGER L.D.M. 9/11/08*,//) 



TAB /M 
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INPUT Of PARAMETERS 



WRITE ( IUMIT2, 1 ) 

FORMAT (//,* 1) DEGREE Or THE CHARACTER 1 ST I C POLYNOMIAL 
C (1<P<13)? ',$) 

READ (IUNIT3,2)P 
FORMAT ( I ) 
IT ( P . LT. 2 )T!f EN 
WRITE (IUNIT2,3) 

FORMAT {' DECREE TOO SMALL I ' ) 
GO TO 10 
END IF. 

IF f P.GT.12)THEN 

WRITE { IUNIT2,4) 

FORMAT (' DEGREE TOO BIG! ' ) 

GO TO 10 

END IF 

WRITE ( IUNIT2 ,201 )P 

FORMAT {* DEGREE OF CHARACTERISTIC POLYNOMIAL: ',121 

CALL PNUL (D£C,P2) I CREATES THE CHAR ACT • POLYN. WITH 

P2(P+l)-l I LAST AND 

P2{1)-1 t FIRST COEFFICIENT!; (ALWAYS EQUAL TO 

WRITE (IUNIT2,5) I THEN , ASKS FOR OTHER COEFFICIENTS: 

FORMAT (/.' INPUT OF THE CHARACTERISTIC POLYNOMIAL.',//. 

C ' POLYNOMIALSS ARE WRITTEN AS:',/, 

C * X0 ¥ Al XI * A2 X2 ♦ ♦ AP-l XP-i *XP',/, 

C * PLEASE GIVE THE RANK OF COEFFICIENTS Al TO AP-l EQUAL TO 1, 
CONE AFTER THE OTHER.',/, 

C ' (THOSE OF DEGREE 0 AND P ARE EQUAL TO I ALREADY)',/, 
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C ' INPUT 0 TO IND I CATC THE EMU OF" THE OPERA I* I ON . ' , / ) 
20 CONTINUE 

WRITE < IUNIT2, 7) 
70 7 FORMAT (* RANK Of A COtrflCIENT EQUAL TO 1 7 ',$) 

READ ( EUNITj,8) II 
8 FORMAT ( I ) 

IT (II.EQ.O) GO TO 25 I MEANS THE END Or THE INPUT 

P2(IJ+1)-1 1A0 IS IN P2(l), Al IN P2(2)... 

CO TO 20 I NEXT COEF. TO INPUT 

25 CONTINUE I INPUT TERMINATED 

C *••**•••*•••*•*•••• 

75 C THE CHARACTERISTIC POLYNOMIAL HAS BEEN GIVEN. VERIFICATION: 

C 

WRITE ( IUNIT2,9500) 
9S00 FORMAT (/) 

CALL ECRIPOL < IUNIT2 , * CHARACTERISTIC POLYNOMIAL i ', DEG, P2 , 25 , 0 ) 
C *• WRITES THE POLYNOMIAL ** 

WRITE (IUNIT2.9) 

20 9 FORMAT (/, ' OR7 :( RETURN- YES , IF NOT, INPUT: N > ',$) 

REP- * ' 

READ ( IUNIT3, 12JREP 

12 FORMAT (A) 

IT ( REP< lit). CQ. ' N* )GO TO 10 I TO STARTING POINT 
c ***•*•***•*• THE WISHED COEFFICIENTS WERE ENTERED. IS THE POLYNOMIAL CORRECT 7 
c ********* WE VERIFY THAT POLYNOMIAL P2 GENERATES A MAXIMUM P.N. 3. 
25 CALL SEQUENCE ( IUNIT2 , DEC, P , P2 , A , SEQl , IND1C) ICOMPUTES THE P.N.S. 

IP ( INDIC.EQ. 1> THEN 

WRITE (6,19) « 
19 FORMAT </,' THE P.N.S. IS NOT MAXIMUM III') 

GO TO 10 I TO STARTING POINT 

END IF 

c ******** TtlE CHARACTERISTIC POLYNOMIAL IS GOOD 
9000 CONTINUE 
30 WRITE { IUNIT2, 13) 

13 FORMAT (/,' 2) NUMBER OF SIMULTANEOUS BITS TO RE ISSUED 
C (N>P-l)7 ',$) 

READ < IUNIT3, 14)N 

IP (N.LT.P) GO TO 30 IN MUST OE AT LEAST EQUAL TO P 

14 FORMAT ( I ) 
WRITE (I UN IT 2, 9) IOK7 
REP-' * 

READ ( I UNIT J , 12 ) REP 
IF ( REP(L:l).EQ.'N')GO TO 30 
c ********••***•**•*•*«****•**»*••*••***• 

C END OF PARAMETERS INPUT. 

c *«»*»*•*»**•****•****••* ********* 

c ********* ************* *••*•*****••**«•• 

40 C 



30 



35 



50 



55 



WRITE ( IUNIT2,9S00) 
WRITE ( IUNIT2,20J) 
20 3 FORMAT tx,'**V"********* A ****** AA ****** 4 *'> 

WRITE (IUNIT4,205* P 

WRITE < IUNIT2,205) P m % 

205 FORMAT <' DEGREE OF CHARACTERISTIC POLYNOMIAL : ',12,/) 

45 CALL ECRIPOL ( IUNIT4 9 ' CHARACTERISTIC POLYNOMIAL ' 

C, DEG, P2, 25,0) 

CALL ECRIPOL (IUNIT2,' CHARACTERISTIC POLYNOMIAL' 
C,DEG,P2,25,0) 

WRITE UUNIT4,17)N 
WRITE (IUNIT2 f 17)N 
17 FORMAT {/.' NB OP SIMULTANEOUS BITS: *,l3,/) 



C * * * * PARAMETERS ARE WOW SUMMARISED ON THE SCREEN AND PUT IN THE TILE. 



c •••• 

C*» *.*••***•*•*******•• 

C **» END OF INITIALISATIONS 



it '3 
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c 

C *mm 

207 
C •* 
2000 



C * 
C * 
C • 

cz 

2001 
CZ 

2003 
C •* 



CZ 

2005 
c •* 
c •* 



WRITE (IUNIT2,207I 

FORMAT (/,X,' COMPUTING ') 

WRITE <IUNIT2,203) 

LOOKING FOR THE SOLUTION • 

KK-l I FIRST LINE AT THE BOTTOM OF THE MATRIX 

CONTINUE 

IF (KK.GT.N)GO TO 2900 I THE END 
IDEP-aK 

NBC-2 !N8 OE COEFF DE ZIC RK AUTRES QUE Z0 
» WE LOOK FOR A POLYNOMIAL PI MULTIPLE Of P2 AND HAVING ONLY 
' NBC-TWO COEFFICIENTS BETWEEN KK AND KKtN-l.SUCII A POLYNOMIAL IS A 
• GENERATOR OF THE P.N.S. SEQL , ALREADY COMPUTED . 

WRITE ( I UNIT! , 2001 )KK,N+KK-l, NBC 

FORMAT (' KK- ',13,' N+KR-l- ',13,' NBC* ',13) 
CALL POL Y AN COE F X ( DEG, Pi , KK,N+KK-l f NBC, P2 , I ND I CP, SEQL ) 

IF ( INDICP. EQ.1)WRITE ( IUNIT1 , 200 3 ) INDICP i FOR TEST OR DEBUGGING 

FORMAT (' INDICP- ',13) 
.* *• ••••••*••••••••••••»••*•••*«» •••*. CAS ou NBC .2 

IF ( INDICP. EQ.l ) THEN I SUCCESS 

Pl(I)-0 I TO LOOK FOR THE OEGREE OF THE 1ST COEF OTHER THAN 

I Z0, ALWAYS EQUAL TO 1 
CALL DEGDE8 (Pi, DEG , I DEGDEB , I ND I C ) I IDEGDEB-DEGREE OF THE FIRST COEFF 
WRITE ( IUNIT1 , 2005 ) I DEGDEB 
FORMAT (* I DEGDEB- ',13) 

pi MAY BE USED UNTIL K- I DEGDEB 
POLYNOMIAL Pi IS COPIED IN THE MATRIX RN WITH THE CORRECT SHIFT 
DO 2100 KKl-KK, I DEGDEB 

I FINISHED 



IF (KKl.GT.N) GO TO 2900 
DO 2200 LL-l.N t Pi COPIED IN RN 
RN (N-KXl+l,LL)-Pl(LL+KKl) 
2200 CONTINUE 

NOM 4— ' RN * 

WRITE ( NOM4 (3:4), '(12. 2)') N-KK 1 1 f I 

CALL ECRIPOLTAB ( I IDEG,NBRN , IUNIT4 , NOM4 , DEG, RN , M-KK 1 »• I , 70 , 0 ) 
2100 CONTINUE 

KK-C DEGDEB* 1 
GO TO 2000 

END IF. I IF INDlCP-0, NOSUCCES5, WE TRY WITH 3, 4, ETC. COEFF 



C 

2300 

CZ 

CZ 



22S0 



C 

2900 
3000 
C »*« 
C *•« 
C *•« 

c 

c **' 



CONTINUE 
NBC-NBC* I 

WniTE ( lUNlTl,200I)KK,NfrKX-l,HBC 

CALL POLYANCOEFX ( DEG r PI , KK , N I . W1C, P2 , I NDl CP . SEQl ) 

IF { INDICP. EQ.l JWRITE ( IUNITl # 200 1 ) INDICP 

IF ( INDICP. EQ.l) THEN tSUCCES 

DO 2250 LL-l,M 

RN ( N-KKtl , LL)-Pl ( LL+KK ) 

CONTINUE I PI PUT IN RN 

NOM4-»RN ' 

WRITE CNOH4I J:4).Ml2.2)')N-KKI*i 

CALL ECRIPOLTAB ( IIDEG.NBRN, UJN1T4 ,NOH4 ,DEG, RN, N-KK l-l , 70 ,0 ). 
KK— KK* I 

CO TO 2000 I NEXT LINE Or THE MATRIX 

END IP t OR NOSUCCESS 

GO TO 2300 I ONE MORE COETF . 

CONTINUE 
CONTINUE 
• * THE END .... 



END OF THE COMPUTATION 



WRITING THE MATRIX IN THE FILE 

WRITE ( IUNIT4, 607 KHAR (12) I FORM FEED 
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IF (N.CT.J2) THEN 

WRITE { IUNIT4, 611) 
611 FORMAT (//,' MATRIX TOO BIG, UNABLE TO WRITE IT... * ,//) 

GO TO 600 

END If 
607 FORMAT (X,A) 

w WRITE ( TUMIT4 , 9500 ) 

WRITE (IUNIT4,207) 

WRITE (IUNIT4,205) P 

CALL ECRIPOL (IUNIT4,* CHARACTERISTIC POLYNOMIAL 
CDEG,P2, 25,0) 

WRITE (IUNIT4,17)N 
WRITE (IUNIT4,203) 
75 WRITE ( IUNIT4,9500) 

WRITE (IUNIT4,601) 
601 FORMAT <//. 10X, 'MATRIX; ' ,/> 

IF (N.LE.24)KE-3 
If (N.GT.24)KE-2 

WRITE ( IUNIT4,605)( { JJ- I )/10 , JJ- i , N) 
605 FORMAT < 5X # <N>( KKE-1> ,X)) 

20 WRITE ( IUNIT4,603)< JJ-1-10M ( JJ-D/10) , JJ-1,N> 

603 FORMAT <5X f <N>< KKE> ),/) 

DO 600 II-1,N 

LIGNE- ' ' 

OO 610 JJ-1,N 

JJl-KE*JJ 

If (RN(It,JJ).EQ.O)LlGNE ( JJl : J Jl,)-' - 9 
25 ir ( RN( X X , JJ | • EQ . X ) LIGNE ( JJl : JJl )«•**•• 

610 CONTINUE 

WRITE (IUNIT4, 617)11-1, LIGNE 
617 FORMAT (X, I2,2X,A75) 

600 CONTINUE 

WRITE ( IUNIT4 ,9500) 
C ********* THE MATRIX IS WRITTEN 
30 C 

C *************** VERIFICATION **************************** 

IP£RMA*X-2**P-1 
NOMDREDEO ITS — 2 • I PERMAX ♦ 1 00 
c • «,****»***•*** VERIFICATION OF THE PARALLEL GENERATOR: ******* 
C *** WE USE POLYNOMIAL 'A' FOR STATE M AND FOR STATE M*l 

C WE LOAD THE GENERATOR WITH THE FIRST N VALUES OF SEQl 

35 DO 5000. II-l.N 

5000 A(N-II + l)-SEQl( II) I TAKE CARE OF THE TIME INCREASING DIRECTION! 

C »•* WE NOW COMPUTE THE VALUES TO PUT INTO B: 

C 

DO 5500 M— 1 , NOMDREDER IT5/N 
OO 5100 tl-UN 
4 0 B(tI)-0 

DO 5200 JJ-1,N 

IF (RN(IZ,JJ|.EQ.t) B(Il)-B(tI).XOR.A(JJ) 
5200 CONTINUE 
5100 CONTINUE 

C ** WE PUT THE VALUES OF 8 IN 5EQ2: 
DO 5300 Il-l,N 
45 SEQ2C (M-l)*N*II)-A<N-ll*l) 

5 300 CONTINUE 

C WE PUT THE VALUES OF B IN 'A* FOR THE FOLLOWING STATE J 

DO 5400 Il-L,N 

5400 Mtn-Bitn 

5500 CONTINUE 

C ************* END OF THE COMPUTATION OF SE02, OBTAINED WITH THE PARALLEL 
50 C GENERATOR, AND COMPARISION WITH SEOl, OBTAINED WITH A SERIES GENERATOR, 

WRITE ( IUNITl,4?03MSEOUKK),KK-l,72) 
4703 FORMAT (' SEOl : ' * X , 72 1 U 

WRITE < IUNITl,570 3 MSEO2{KK) ,KX-1,72) 
5703 FORMAT (' SEQ2 i ' ,X, 7211 ) 



55 
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CO 6000 r I-l, 2MPERMAX 
5 IP (S£QU H).NE.SEQ2(n)) GO TO 6200 

6000 CONTINUE 

WRITE < IUNIT4,6001) 

6001 FORMAT (//,' VEBtriCATlON O.K. Ill') 
GO TO 10000 

6200 WRITE (IUNIT4,6201) 

6201 FORMAT (' THE P.N.S. ARE DIFFERENT. Ill SOMETHING WRONG.') 
10 100O0 CONTINUE 

WRITE (IUNIT2,209) 
209 FORMAT (//, ' JOB TERMINATED . RESULTS IN GPSA.DAT' ,//) 

END 

SUBROUTINE SEQUENCE ( IUN IT , DEC , P , P2 , A. , SEQ 1.1 ND I C > 
C *• ON GENERE LA SEQUENCE QUE FOURNIRAIT LE SHIFT REGISTER A 
75 C DE P BASCULES ET DE POLYNOMS CARACTER ISTIQUE P2 . 

C DEUX PERIODES DE LA SEQUENCE KAXIMALE SONT RANG EES DANS SEQl. 

C SI LA PERIODS N' EST PAS MAJCIHALE, L ' I ND I CAT EUR EST HIS A I 

C LES MESSAGES SONT ECRITS SUR IUNIT 

C 

C ** GENERATES THE SEQUENCE SEQl FURNISHED BY THE SHIFT REGISTER A, 
20 C »» Of P STAGES AND CHARACTERISTIC POLYNOMIAL P2 

C *• TWO PERIODS ARE COMPUTED . tF THE P.N.S. IS NOT MAXIMUM, INDIC-l 
C ** MESSAGES ARE WRITTEN ON IUNIT1 
INTEGER DEG,P 

INTEGER P2(l) # A(l) ,SEQl(l) 
INDIC-0 

C LOADING 'A» WITH THE SEED 

25 CALL PNUL( DEG, A) 

DO 4100 II-2,P»-1 
410O A(II)-l I DONE 

C ********* ALGORITHM ************* 
IPERMAX-2**P-1 I PERIOD 

NOMBREDEBITS-2*IPERMAX+100 I NOMB REDED ITS > 2 * I PERMAX 

DO 4S00 MM— I , NOMQREDEBITS I MORE THAU TWO PERIODS. 

30 C *** ON CALCULE L ' ELEMENT QUI VA ENTRER DANS LE REG I STRE . C • EST AC I ) 
A<1)-AIP*1) I ALWAYS CONNECTED 

DO 4200 II-2,P I IT AMI) CONNECTED. . . 
IF (P2fII).EQ.l) Afl)-A(l) .XOR.A( II) 
4 200 CONTINUE 

C ** COLLECTING THE ISSUED OCT AND SHIFTING 
SEQMMM)-A(Ptl) I COLLECTING 
35 DO 4300' II-P+l,2,-l I SHIFTING 

4300 At II )-A< II-l) 
4S0O CONTINUE 

c ******* ON VERIFIE LA PERIODS • PERIOD VERIFICATION. 

DO 4600 I PER* 1,1 PERMAX* 1 

DO 46S0 II-l,IPER+20 
40 IF <SEQI< Il+IPER) .NE.SEQK II) ) GOTO 4600 

4650 CONTINUE 

GO TO 4700 
460O CONTINUE 

4700 CONTINUE 

WRITE (IUNIT, 470t)lPER,IPERMAX 

4701 FORMAT <//,' PERIOD OF SEQl: 
45 C 16, » MAX PERIOD : ',16) 

WRITE ( IUNIT, 4703 ) ( SEQl ( KK ) , KK- 1,72) 
4703 FORMAT (' SEQl :',X,72Il) 

I F ( I PER . NE . I PERMAX ) 1ND I C- 1 

RETURN 

END 

50 C • 

C M'.M ***** ****************************** 

c ******** POLY 2JDP.FOR 

c *********** **T***** ******************** ******** ******** * ************ 

C •* BI8LIOTIIEQUE D'OPERATIONS SUR LES POLYNOME5 DOMT LSS COEFFICIENTS SONT 
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c •* oes cmtiers mooulo 2 

c 

C L I 5TE DCS SOUSHOUTINES DC TRAITEMENT DE POLYNOMES A COEFFICIENTS 
° C SUB LE CORPS i 0.1) 
C 

C LES POLYNOMES SONT ORDONNES PAR DEGR E CROISSANT 

C 

C COPIE <DEG,Pl,P2T :RECOPIE Pi DANS P2 

C DEGRE { POL, DEG, DEGM, INDTC) : DONNE LE DEC RE DE POL 

C DEGDEB ( POL, DEG, DEGB , 1NDIC : DONNE LE PREMIER TERME NON NUL 

75 c PNUL <OEG,P) : CREE LE POLYNOME P I DENT I QUCftENT NUL 

C PUNIT { DEG , P ) : CREE LE POLYNOME EGAL A 1 

C ECRIPOL (IUNIT,NOM,OEC,P,DEGMAX, INDIC) 

C :ECRIT LE POLYNOME P SUR L' UNITE IUNIT, 

C PRECEDE D'UN TITRE *NOM'. JUSO'A DEGMAX 

C SI INDIC-O, LES TERMES SONT GROUPES 

C SI INDIC-l, LA PLACE DCS TERMES NULS 

2Q c EST REMPLACEE PAR DES 8LANCS . 

C ECR I POLTAB ( I IDEG, NBRN, lUNtT, NOM, DEG, P , I . DEGMAX, INDIC) 

C : COMME ECRIPOL, MAIS POUR UN TABLEAU 

C DE POLYNOMES, OONT ON ECR IT LA LIGNE I. 

C POLYANCOEFX i CHERCHE LES MULTIPLES D'UN POLYNOME 

C : QUI ONT UN NOM8RE DONNE DE COEFFICIENTS 

£**••*•»*•*•*»**••»**•*•*»•*****••**•*****•••***•*»*•****• 
c REMARQUES GENERALES . REMARKS. 

25 C LES POLYNOMES SONT TOUS ECRITS DANS UN VECTEUR A DEG POSITIONS 
C EN PART ANT DU DEGRE NUL. 

C POLYNOMIALS ARE WRITTEN IN AN ARRAY, STARTING TRON 0 DEGREE TERM . 

C LE DEGRE MAXIMAL TRAITABLE EST DONC DEG- 1 . ATTENTION AUX DEBOROEMENTS I M 

C LES POLYNOMES AUXILIAIRES DONT LES DEGRCS NE SONT PAS PASSES EN ARGUMENT 

C SONT DIMENS IONNES A 256. 

C 

30 C 

SUBROUTINE DEGRE ( POl*, DEGN,DEG, INDIC) 
C DONNE LE DEGRE DU POLYNOME POL ECR IT DANS DEGM CELLULES 
C *• DEGM-DEGRE MAX* 1 , DEG-DEGRE DU POLYNOME 

C *• LOOKS TOR THE DEGREE Or THE POLYN . INDIC- 1 IT THE POLYN. IS NULL 
INTEGER DEGM , DEG , POM M 
35 INDIC-0 

DO 200 II-DECM,l,-l 
IT (POMU).E0.1)GO TO 210 
200 CONTINUE 

IF (POL(l).EQ.O) THEM 

DEG— 1 UNDICATION DE DEGRE NUL 

C WRITE (6,1) 

40 i FORMAT ( X, ' DEGRE I POLYNOME IDENTIQUEMEMT NUL I * ) 

INDIC- I 
C GO TO 10000 

END ir 
210 OEG-lI-i 
10000 CONTINUE 

RETURN 

45 ENO 
C 

SUBROUTINE DEGDEB < POU, DEGM, DEGB , INDIC) 
C DEGM-DEGREMAX , DEGB-DEGRE DU 1ER TERME OU POLYNOME 

C INDIC- I SI POLYNOME NUL 

C LOOKS TOR THE FIRST TERM Or THE POLYN. 

5 0 INTEGER DEGM, DEGB ,POL(l) 

INOIC-O 

DO 200 I t-l , DEGM 
IF (POL( II) .EQ.DGO TO 210 
200 CONTINUE 

IF ( POL ( DEGM ) • EQ . 0 ) THEN 
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n-7 
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c 


WFITE (6, i) 


i 


ronMAT ( X , ' DEGDEB : POLYNOME ILrEMTtQULMENT NOLI') 




It 0IC-1 


c 


GO TO 10000 




END if 


210 


OEGB-II-1 


10000 


CONTINUE 




RETURN 




END 



c ******* 

c *•••**•**••*****•***•*«•**.«*..*. 

SUBROUTINE COPIEX (DEG,P1,P2) 
C ** RECOPIE PI DANS P2 TERME A TERME. COPIES Pi INTO P2 

INTEGER DEG 

INTEGER PU i) ,P2( 1) 

DO 100 I t - I , DEG 

P2(II)-P1(II) 
100 CONTINUE 

RETURN 

END 

c ************ ••*•«•************»«*•** ******* 

SUBROUTINE PNUL ( DEG , P ) 
C ** CREE LE POLYNOHE NUL OE DEGRE MAX DEG • CREATES THE NULL POLYNOMIAL 

INTEGER DEG 

INTEGER PU) 

DO 100 11*1 , DEG 
100 F( II )»Q < 

RETURN 

END 

c *»••**** ************ ***************** 

SUBROUTINE PUNIT <DEG,P) 
C **••♦• CREE LE POLYNOHE OB DEGRE 0. CREATES THE CONSTANT POLYN. •**•••* 

INTEGER DEG 

INTEGER PU) 

CALL PNUL (DEG,P) 

P( 1)-1 

RETURN 

END 

c A***************************************** 

C ********************************************************* 

SUBROUTINE ECRIPOL ( IUNTT, NOH , DEG # P # DEGMAX , INDIC) 
C •* SI INDIC-0,TERMES BLOQUES.Sl INDIC-1 , ESPACES OE 4 BLANCS POUR CIIAQUE 

C ** TERME NUL 

C ** IUNIT:L'UNITE LOG I QUE SUR LAQUELLE ON ECRIT. 

C •* WRITES A POLYNOMIAL P. IT INOIC-O, NULL TERMS ARE DC SCARDED. I f INOlC-l, 
C •* FOUR BLANKS ARE LEPT TOR EACH BLANK TERM. DEGMAX LIMITS THE NUMBER Or 
C ** WRITTEN TERMS. NOM IS THE NAME OF THE POLYNOMIAL, WHICH MAY BE WRITTEN. 

CHARACTER* ( *)NOM 

CHARACTER* 80 LIGNE 

INTEGER DEG t DEGMAX , DEGMAX 1 

INTEGER PU) 

CALL DCGRE ( P , DEG, DEGMAX 1 , INDIC2 ) 

L-LEN (NOM) 

LIGNE-' 9 

LIGNE UtD-NOM 

LIGNE (LUsLO)-' : • 

IDEF-L+4 

c ************************************************ 
DO 100 1 1*1 .DEGMAX 

IF (PUI). EQ.0. AND. INDIC. EQ.l) GO TO 120 ION SAUTE 4 BLANCS 

IT (P(n).EO.U THEN 11 
LIGNEC IDEP: IDEP)-'Z' 

IF ( (I l-l) .LE.9)WRITE( LIGNE( I DEPU : IDEFU ) 9 • (tl )' )! I- I 
IF ( ( I I - I ) . GT • 9 . AND , ( I l-l) . LT . 100) 
C WRITE( L IGNE{ IDEPf 1 : IDEPf 2 ) , '( I 2 ) • >II-l 
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IT t(-II-i..GT. 99. ANlKtll-ll-LT. 1000) 
C WRlTEt LlGNE< I DEP * 111 DEP 3 ) ,'(I3)')II-1 
120 IDEP-IDEP+4 

lf((lDEP^8).GT.eO) THEN 12 

LIGNEt IDEPi IDEP+3)-' ' 

WRITE ( tUNXT,i0l> 
iOt FORMAT (X,* RESULT 13 CUT') 

GO TO 1000 

END IP J* 
END If 11 

100 CONTINUE 
1000 CONTINUE 

£♦»•*•••••••••*•**••••••*•**•**• *'* 

WRITE ( IUNIT, 103 )LIGNE 

101 FORMAT (X,A80) 
RETURN 

C ECRIT LE POLYNOME DC RANG I OU TABLEAU DE POLYNOMES P 

5 Itl.i^ilUtSri'fSlWWM. Or RANK I. TAKEN IN A TABLE OP POLYNOMIALS 
CHARACTERMMNON 
CHARACTER* B0 LIGNE 

INTEGER OEG * DEGMAX , DEGMAX I 

INTEGER P{ NBRN r 1 1 DEG ) # A ( 300) 

DO 10 tl-l,DEG 

10 CALL^ ECRIPOL^ < I UN IT, NOM , DEG , A , DEGMAX , INOIC) 

RETURN 

% EXAMPLE, ir N-J , U*"*£ - ' 0 «.«. ^ b. found. 

ESSS SfiWBTsXSKs" 

SSSSr !«lt« JR 1 ** 1 ."SITIOMS Or THE INDICES. 
r»LL MGH ( P. OEG. IOECP. IMDICIU 
M^X-2"IOECP-l IPEMOD Of THE P.M. I. 

" - , , u „ , r _i ni CAS DC succes. , 

IHOtC-0 „^nea-» THE TIHST POStTtOM WILL DE atl „_„„„..,- 

1DEP-0EGDEB ! *nf u TMKS ARt POS I T IONNED BLOCKED TOGETHER. STARTING 

BO 10 II-l.N | f ROM DEG0EB . THEY HI LL BE MOVED APTER EACH TEST. 
C IHOlCeHD-tOEPMt* 0 ? ?H01CEUI) IS THE POSITION Of TERM II 

5«5S««i7W-i i^.r^o^r'" TERMS * T,, " R PLACES 

120 CONTINUE TEST 

% z "« t„e .■ * or s E0 . 



31 



EP 0 397 079 A2 



20 



00 2 100 [ECH-L*OEG . , UMAX ► DtCMAX t U.LH I lie *o\nfv uf owe oil 

UESV-SEQ (I ECU J I OF Til E c.N.S. SEQ 

DO 2000 ir-i,N 

ITEST - ITEST. XOR. SEQ< IECH-INDICEI tl)«-l> 
10 2000 CONTINUE 

IF UTEST.NE.0) GO TO 2200 t INTERRUPTS THE TEST AS SOON AS A 

C I DISCREPANCY APPEARS * 

210 0 CONTINUE 

GO TO iOOO I SUCCESS I ! I 

C •* SINON, ON ESSAIE DE DECALER LES INDICES 
C IF NO SUCCESS, THE TERMS OF THE POLYN , ARE SHIFTED 

75 2200 CONTINUE 

c END Or TEST »*••**•*•••••♦•*•*••••••*•* 

DO 200 II-i.N-l 

ir (INDICE(Il) .LT.INDICE<II*l)-l) THEN t THERE IS ROOK TO SIIITT THE TERM 
INOICE <II)-INDICE{II)*1 tIT IS DONE 

DO 230 LL-1,II-1 I PRECEEDINQ TERMS ARE PLACED AGAIN AT THE 

230 IND ICE( LD-LL+IDEP I BEGINNING 

GO TO 100 I TO THE TEST 

ENO IF 

200 CONTINUE I WHEN HERE, ALL TERMS ARE BLOCKED AGAINST THE 

C I NTH TERM, WHICH MUST BE SHITTED 

INDICE(Nl-INDICE(N)+l I TERM N IS SHIFTED 
C FIN 

IF ( INDICE( N) .GT .DEGMAJC+ 1 ) THEN ISTOP'.N GREATER THAN DEGMAX 
INDIC-O INO SUCCESS AT ALL. EXIT 

25 GO TO 10000 

ENO IF 

C *»4******4*»» 

DO 210 LL-l,N-l t OTHER TERMS ARE PLACED AGAIN AT THE BEGINNING 
INDICEf LL ) • IDEP+LL 
210 CONTINUE 

GO TO 100 ITO THE RETURH POINT 

30 1000 CONTINUE ISUCCESSIM! 

INDIC-1 
10000 CONTINUE 
RETURN 
ENO 

35 
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SAMPLE TERMINAL LlSTiMfi 
-RUNNING PGM GPSA - E 



TABLE 12- 



$ RUN CPSA_E 

PARALLEL PSEUDONOISE SEQUENCIES GENERATOR. 

COMPUTATION OF THE TRANSITION MATRIX. ( PROGRAM CPSA) 
CEORGES ROGER L.D.M. 9/11/88 



25 1) DECREE OF THE CHARACTERISTIC POLYNOMIAL (1<P<13>? 7 

DEORE 0U POLY NOME CARACTER 1ST IQUE : 7 

INPUT OF THE CHARACTERISTIC POLYNOMIAL. 

POLYNOMIALSS ARE WRITTEN AS: 

XO ♦ Al XI ♦ A2 X2 ♦ ♦ AP-l XP-1,*XP 

30 PLEASE GIVE THE RANK OF COEFFICIENTS Al TO AP- 1 EQUAL TO I, ONE AFTER THE OTHE 

R. 

(THOSE OF DEGREE 0 AND P AnE EQUAL TO 1 ALREADY ) 
INPUT 0 TO INDICATE THE END OF THE OPERATION. 

RANK OF A COEFFICIENT EQUAL TO 1 ? 6. 
flJVNK OF A COEFFICIENT EQUAL TO 1 ? O 

35 

CHARACTERISTIC POLYNOMIAL; : ZO 26 7.7 
K? : (RETURN«YES, IF NOT, INPUT: N ) Y 

PERIODE TROUVEE POUn SEQ L : 127 PERIOOE MAX: 127 

40 EQl: 1 11 1 1 1 1OOO0OO 100000 1 1 0000 101000 111 1001000 101 100 1 I 10 1 0 100 1 1 it 1 0 10000 t I 100 

) NUMBER OF SIMULTANEOUS B ITS (H>P-l)? ^ * 

Kl : (RETURN* YES , IF MOT, INPUT: H ) Y^ 

DECREE OF CHARACTERISTIC POLYNOMIAL : 7 
CHARACTERISTIC POLYNOMIAL : 20 Z6 27 
3 OF SIMULTANEOUS BITS: 8 



50 

COMPUTING 

soi"iiitiiioooaooioooooi looooioiooonuooioooioi 1001 tioioioouu 10100001 uoo 
SQ2 : . in in tooooooioooooiiooooioioooiiuooioooioiiooiiioio 1001 111 10 10OO01 1100 

55 JOB TERM I HATED. RESULTS IH OPSA.DAT 

TAB*12-i 
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$ TY OPSA.DAT 

PARALLEL PSEUDONOISE SEQUENCIES GENERATOR. 

COMPUTATION OF THE TRANSITION MATRIX . ( PROGRAM CPSA ) 
GEORGES ROGER L . D . M . 9/11/88 



10 



75 



20 



DEGREE OF CHARACTERISTIC POLYNOMIAL : 7 

CHARACTERISTIC POLYNOMIAL : 20 26 Z7 
NB OF SIMULTANEOUS BITS: 0 



RN07 : 


25 


26 


RN06 : 


24 


25 


RN05 : 


23 


24 


RN04 : 


22 


23 


RN03 : 


: Zi 


Z2 


RN02 . 


20 


21 


UNO I : 


25 


27 


RNOO : 


24 


26 



25 



30 



35 



pECREE OF CHARACTERISTIC POLYNOMIAL : 7 
CHARACTERISTIC POLYNOMIAL : 20 26 27 

HB OF SIMULTANEOUS BITS: B 
I********************************* 



40 



45 



50 



55 



MATRIX: 

00000000 
0 1 2 3 4 5 6 7 

0 - - + 

1 + - + 

2 f*--_--._ 

3 - + *- -- _- 

4 -. + + .- _- 
5 



TAB. 12-2 
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TABLE 13 > 

DEGREE Or CHARACTERISTIC POLYNOMIAL. 

CHARACTERISTIC POLYNOMIAL : 10 16 
NB Or SIMULTANEOUS BITS: 6 



MATRIX: 

00000000 
01234567 

0 ---_ + - + - 

1 _•--_ + - + 

2 + + 

3 - + + -.--- 

4 -- + + -»-- 

5 - 

S • -•. + + .- 

7 - 



VERIFICATION O.K. Ill 



TAB. 13-1 
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PARALLEL PSEUDONOISE SEQUENCIES GENERATOR . 

COMPUTATION Of THE TRANSITION MATRIX. ( PROGRAM OPSA) 
PEORviES ROGER L.D.M. 9/11/8 8 



DEGREE OF CHARACTERISTIC POLYNOMIAL : 7 

CHARACTERISTIC POLYNOMIAL : ZO Z6 Z7 
NB Or SIMULTANEOUS BITS: : B 



20 



RN07 \ 


ZS 


Z6 


RN06 J 


24 


ZS 


RN05 ! 


Z3 


Z4 


RN04 ' 


22 


Z3 


RN0 3 


: Zl 


Z2 


RN02 


: ZO 


Zi 


RN01 


x Z5 


Z7 


RMOO 


: Z4 


Z6 



25 



30 



35 



40 



45 



TAB 13-2 

50 
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DECREE Or CHARACTERISTIC POLYNOMIAL, : 7 
CHARACTERISTIC POLYNOMIAL : 20 26 27 
NB Or SIMULTANEOUS BITS* (fT> 



MATRIX: 

0000000000111111 
0123456789012345 

0 _-•---_•_ + _- - + _- . 

1 - -- -- -- - - 

2 - - + 

3 « + 

4 + 

5 + 

7 - -- -- -- -- -- 



6 ---- + - + 
9 + _ 

10 ♦ 

11 - + + 

12 ._• + + --- 
13 

14 _-_- + + -. 
15 



VERiriCATION O.K. tit 



t 

TAB. 13-3 
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PARALLEL PSEUDONOISE SEQUENCIES GENERATOR 

COMPUTATION Of THE TRANSITION MATRIX. 

GEORGES ROQER L . D . M . 9/11/B8 



DEGREE Or CHARACTERISTIC POLYNOMIAL : 7 
CHARACTERISTIC POLYNOMIAL : Z0 Z6 Z7 



NB Or SIMULTANEOUS BITS: 




RN15 : 


25 


26 


RN14 : 


24 


25 


RN13 : 


23 


24 


HN12 5 


22 


23 


RN11 ! 


21 


22 


RN10 s 


20 


Zl 


RN09 : 


25 


27 


RN06 i 


24 


26 


RN07 


l 23 


25 


RN06 


: 22 


24 


RN05 


: 21 


23 


RN04 


: 20 


22 


RN03 


: 211 


Z15 


RN0 2 


: 210 


Z14 


RN01 


: 29 


Z13 


RN00 


: Z0 


Z12 



TAB. 13-4 



38 



EP 0 397 079 A2 



PARALLEL PSEUOONOISE SEQUENCXBS GENERATOR, 

COMPUTATION Or THE TRANSITION MATRIX, { PROGRAM CPSA) 
GEORGES ROGER L.D.M. 9/11/99 



DEGREE Or CHARACTERISTIC POLYNOMIAL I 7 
.JUlAJU^KRJJLTJJ^.J&ftrywoRtAfc..^ 20 Z6 27 
NB OP SIMULTANEOUS BITS? 24 



RN23 


l Z5 


26 


RN22 


1 14 


25 


RN21 


1 23 


24 


RN20 


1 It 


23 


RNl 9 


\ 21 


£3 


RNie 


I 20 


21 


RN17 


1 25 


27 


RNl 6 


1 Z4 


26 


RN1S ! 


1 23 


z£ 


RN14 


I 22 


24 


RNl 3 i 


I 21 


23 


RN12 


I 20 


Z2 


RN11 


i 211 


215 


RNiO i 


I 210 


214 


RH09 i 


i 29 


213 


RN05 I 


t 28 


212 


RN07 i 


i 27 


211 


RN06 i 


I 26 


210 


RHOS I 


i Z5 


29 


UNO 4 I 


i 24 


2* 


RN03 1 


I 23 


27 


RN02 I 


l 22 


26 


RN01 i 


1 Zl 


25 


RNQD i 


I 20 


24 



TAB. 13-5 
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10 



75 



20 



25 



45 



DEGREE Or CHARACTERISTIC POLYNOMIAL » 7 
CHARACTERISTIC POLYNOMIAL 1 ZO Z* 17 
NB OP 6 I MULTANEOUS BITSl 24 

********************************** 



50 



55 



MATRIX I 



4 

•> - - - 

30 6 - - 

7 - - - 

0 - - - 

9 - - - 

10 - - - 

11 - - - 

12 * - ♦ 

35 13 - ♦ " ♦ 



14 + - + 

15 
16 
17 

10 + + 



19 - + ♦ - - 

40 20 - - + ♦ " 
21 - - - 4 + 



22 
23 



♦ - - - 



VERirXCATION O.K. Ill 



TAB 13-6 
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PARALLEL P8EUDONOI3E SEQUENCIES GENERATOR, 

COMPUTATION Or THE TRANSITION MATRIX. ( PROGRAM GPSA) 
GEORGES ROGER L.D.rt. 9/11/BS 



DEGREE OF CHARACTERISTIC POLYNOMIAL i 7 

CHARACTERISTIC POLYNOMIAL t 20 Z$ 11 
Nfi OP SIMULTANEOUS BITSj 32 



RN31 


c 


25 


26 


RN3 0 


t 


24 


2b 


RN29 


t 


23 


24 


RN28 


t 


Z2 


23 


RN27 


t 


Zl 


22 


RN26 


t 


20 


Zl 


RH25 


I 


25 


27 


RN24 


1 


24 


26 


RN23 


• 
• 


23 


25 




r 


22* 




RH21 


5 


21 


•p s 
2 J 


nri * v 






4b 4b 


RN19 




Zll 


215 


RNie 




210 


214 


RH17 




2* 


Z13 






28 


112 


RN1S 




Z7 


Zll 


RN14 




26 


Z10 


RW13 




ZS 


29 


RN12 




Z4 


20 


RN11 




23 


Z7 


RNlO 




Z2 


26 


RN09 




Zl 


26 


UNO 8 




20 


24 


RN07 




Zl 


216 


RN06 




20 


217 


RN05 




Zfl 


226 


RN04 




Z7 


225 


UNO 3 




26 


224 


RN02 




ZS 


«J 


RN01 




24 


Z22 


KNOO 




Z3 


221 



TAB 13-7 
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DEGREE Or CHARACTERISTIC POLYNOMIAL I 7 
CHARACTERISTIC t POL YHOMIAL I 20 26 27 
NB Or SIMULTANEOUS BITS I 32 



MATRIX 1 

00000000001111111111222222222233 
01234667890123456789012346678901 

0 ♦ - + 

1 4 + 

2 - - - + - 

3 - * ♦ 



5 

6 + 



7 - + 
B + - 



9 - + - 
10 * 



11 + 

12 



13 

14 

15 — - - 

16 

17 

10 

19 

20 + - + 

21 - * - 
22 ♦ 

23 

24 

25 - - - 
; < - 



27 - ♦ + 

20 ++ 

29 ♦ ♦ - - 

30' ♦ ♦ - 

31 ♦ + 



VERITICATIOH O.K. II I 



TAB. 13-8 
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PARALLEL PSEUDONOISE BEQUENCIES QENERATOR. 

COMPUTATION OP THE TRANSITION MATRIX. { PROGRAM GPSA) 
GEORGES ROGER L.D.M. 9/il/BB 



DEGREE OF CHARACTERISTIC POLYNOMIAL t 7 

CHARACTERISTIC POLYNOMIAL i 20 26 Z7 
NB Of SJMULTANIOUB BITS I 64 * 



RN63 


\ zB 


26 






nN62 


l Z4 


25 






RN61 


! 23 


24 


Kr* 13 1 


Zt 


RN60 


I Z2 


23 


RHB* 


t ~ 1 




RN12 t 


26 


RN58 


! Z0 


21 


RNll I 


27 


RN57 


I Z5 


27 


RNIO I 


26 " 


RN56 


i 24 


26 


RN09 I 


25 


RN55 


I 23 


25 


RN08 : 


24 


RN54 


: 22 


24 


RN07 : 


23 


RN53 


1 zl 


23 


hN06 1 


22 


RN52 
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It will thus be seen that the object set forth above and those made apparent from the preceding 
description are efficiently obtained, and since certain changes may be made in the above construction and 
methodology without departing from the scope of the parallel pseudo-random generator invention, it is 
intended that all matter contained in the above description or shown in the accompanying drawings shall be 
interpreted as illustrative and not in a limiting sense. 

It is also to be understood that the following claims are intended to cover all of the generic and specific 
features of the parallel pseudo-random generator invention herein described, and all statements of the 
scope of the invention which, as a matter of language, might be said to fall therebetween. 



Claims 

1. A parallel pseudo-random generator for emulating a serial pseudo-random generator that generates 
serial outputs such that the next serial output value is based upon an Exclusive OR combination of at least 
two preceding serial output values the maximum preceding serial output value defined as the Pth 
preceeding serial output value, where P is an integer greater than one; comprising: 

A) at least P latches, each latch having an output having a logic value 1 or 0 and an input operable 
upon receipt of a clock signal, for receipt of data for controlling the next logic value on the latch output; 

B) at least P Exclusive OR gates, each having at least two inputs and one output, each Exclusive OR 
gate output connected to a corresponding input of one latch so as to define the next value of the latch 
output upon receipt of the next clock signal; and 

C) means for connecting each input of each Exclusive OR gate to one latch output so that the output 
of each Exclusive OR gate represents the corresponding next value of the latch to which This Exclusive Or 
gate output is connected. 

2. A parallel pseudo-random generator as defined in Claim 1, wherein the serial Exclusive Or 
combination defining the serial pseudo-random generator determines its next output value based upon the 
sixth and seventh preceeding serial output values (P = 7). 

3. A parallel pseudo-random generator as defined in Claim 2. wherein the number of latches is eight 
the latches having corresponding outputs QO through Q7, and wherein the corresponding Exclusive OR 
gates Ex0-Ex7 each having their output connected to the corresponding latch input, have their inputs 
connected to the following latch outputs: 

ExO inputs connected to Q4 and Q6 
Ex1 inputs connected to Q5 and Q7 
Ex2 inputs connected to QO and Q1 
Ex3 inputs connected to Q1 and Q2 
Ex4 inputs connected to Q2 and Q3 
EX5 inputs connected to Q3 and Q4 
Ex6 inputs connected to Q4 and Q5 
Ex7 inputs connected to Q5 and Q6 

4. A parallel pseudo-random generator as defined in Claim 1, wherein the serial Exclusive OR 
combination defining the serial pseudo-random generator combines the sixth and seventh preceeding serial 
output; wherein the number of latches is sixteen, the latches having corresponding outputs QO through Q15, 
and a width of the pseudo-random generator is equal to 16 and further wherein the corresponding sixteen 
Exclusive OR gates ExO - Ex15 each having their output connected to the corresponding latch input, have 
their inputs connected to the following latch outputs: 

ExO inputs connected to Q8 and Q12 
Ex1 inputs connected to Q9 and Q13 
Ex2 inputs connected to Q10 and Q14 
Ex3 inputs connected to Q11 and Q15 
Ex4 inputs connected to QO and Q2 
EX5 inputs connected to Q1 and Q3 
Ex6 inputs connected to Q2 and Q4 
Ex7 inputs connected to Q3 and Q5 
Ex8 inputs connected to Q4 and Q6 
Ex9 inputs connected to Q5 and Q7 
Ex10 inputs connected to QO and Q1 
Ex1 1 inputs connected to Q1 and Q2 
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Ex12 inputs connected to Q2 and Q3 
Ex13 inputs connected to Q3 and Q4 
Ex14 inputs connected to Q4 and Q5 
Ex15 inputs connected to Q5 and Q6. 
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